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Abstract. The finite section method is a classical scheme to approximate the 
solution of an infinite system of linear equations. We present quantitative esti- 
mates for the rate of the convergence of the finite section method on weighted 
€ p -spaces. Our approach uses recent results from the theory of Banach algebras 
of matrices with off-diagonal decay. Furthermore, we demonstrate that Banach 
algebra theory provides a natural framework for deriving a finite section method 
that is applicable to large classes of non-hermitian matrices. An example from 
digital communication illustrates the practical usefulness of the proposed theo- 
retical framework. 

1. Introduction 

Many of the concrete applications of mathematics in science and engineering 
eventually result in a problem involving linear operator equations. This problem 
can be usually represented as a linear system of equations (for instance by dis- 
cretizing an integral equation or because the operator equation is already given on 
some sequence space) of the form 

(1.1) Ax = b, 

where A is an infinite matrix A = {a,ki)k,i& and b belongs to some Banach space 
of sequences. Solving linear equations with infinitely many variables is a problem 
of functional analysis, while solving equations with finitely many variables is one 
of the main themes of linear algebra. Numerical analysis bridges the gap between 
these areas. A fundamental problem of numerical analysis is thus to find a finite- 
dimensional model for (jl.ljl whose solution approximates the solution of the original 
infinite-dimensional problem with any desired accuracy. This problem often leads 
to delicate questions of stability and convergence. 

A simple and useful approach is the finite-section method [11, 17]. Let 

Pnb = (..., 0, 6_ n , 6-n+i, ■ ■ ■ , b n -i, b n , 0, ... ) 

be the orthogonal projection onto a 2n + 1-dimensional subspace. We set 

(1.2) A n = P n AP n and b n = P n b , 
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and try to solve the finite system 



for properly chosen n. The crucial question is then: What is the relation between 
the numerical solution x n and the actual solution x? 

This problem has been analyzed in depth for the case of convolution operators 
and Toeplitz matrices in the pioneering work of Gohberg, e.g. see [11]. Important 
generalizations and extensions in the Toeplitz setting can be found in [4,5]. Rabi- 
novich et al. derive necessary and sufficient conditions for the convergence of the 
finite section method in terms of the so-called limit operator [20], which does not 
necessarily require any Toeplitz structure. These conditions, while intriguing, are 
not always easy to verify in practice. 

A general theory for the approximation by finite-section is based on the power- 
ful methods of C*-algebras and has been developed by Bottcher, Silbermann, and 
coworkers, see for instance [4,17]. Their framework leads to many attractive and 
deep results about the applicability of the finite section method as well as other 
approximation methods. William Arveson goes a step further and concludes that 
"numerical problems involving infinite dimensional operators require a reformula- 
tion in terms of C*- algebras" [1]. However, C*-algebras have some limitations. It 
was already pointed out in [17] that C*-algebra techniques do not yield any infor- 
mation about the speed of convergence of the finite section method. An answer to 
this question is obviously not only of theoretical interest, but it is important for 
real applications. For instance, we want to choose n in ()1.3|) large enough to get a 
sufficiently accurate solution, but on the other hand, n should be small enough to 
bound the computational complexity which in general is of order 0(n 3 ). Theorems 
about the speed of convergence will give a quantitative indication for how increas- 
ing n will impact the accuracy of the solution. Some results about the speed of 
convergence for the special case of Toeplitz matrices can be found in [12,22,25,26]. 
In [12] the convergence in the £ p -norm (1 < p < oo) is analyzed. 

In this paper we present a thorough analysis of the convergence of the finite 
section method for positive definite matrices as well as for non-hermitian ones. 
Specifically, we solve the following problems. 

(a) We study the finite section method on weighted £ p -spaces. If the input vector 
b belongs to a weighted space 0^, then, under suitable assumptions on the matrix 
A, the finite section method converges in the norm of i v m . 

(b) We obtain quantitative estimates for the rate of convergence of x n to x in 
various weighted £ p -norms. 

(c) We define a modified version of finite sections, the non-symmetric finite sec- 
tion method, and show that this method converges also for non-symmetric matrices. 
The finite section method for non-symmetric matrices raises a number of rather 
difficult questions and has motivated a large part of [17]. Even for the classical 
case of Laurent operators (Toeplitz matrices) our approach enlarges considerably 
the class of matrices to which the finite section method can be applied. 

As we work with Banach spaces of sequences, the methods will be taken from 
the theory of _B*-algebras (involutive Banach algebras) instead of C*-algebras which 
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suit only Hilbert spaces. The key property of the matrices A is their off-diagonal 
decay; we will rely heavily on recent results from the theory of Banach algebras of 
matrices. In fact, an important technical part of our analysis is to establish a finite 
section property of infinite-dimensional matrix algebras. 

The paper is organized as follows. In Section 2 we recall the well known proof for 
the convergence of the finite section method for positive invertible matrices and take 
it as a model for more general statements. In Section 3 we introduce several Banach 
algebras of infinite matrices and collect their fundamental properties. Section 4 
is devoted to the notion of inverse-closedness and spectral invariance in Banach 
algebras and their relation to the finite section method. In Section 5 we establish 
the convergence of the finite section method on weighted £ p -spaces, in Section 6 
we derive quantitative estimates. In Section 7 we investigate a version of the finite 
section method for non-symmetric matrices, and in the final Section 8 we briefly 
discuss an application to wireless communications. 

2. Convergence of the finite section method 

It is well known that for positive definite matrices the finite section method works 
in principle, see, e.g., [17]. The proof is instructive and exhibits what is necessary 
for an understanding of the finite section method. 

Recall that if A is an algebra, then the spectrum of an element A G A is defined 
to be the set (74(A) = {A G C : (A — XI) is not invertible}. If the algebra is B(H), 
the bounded operators on some Hilbert space, we usually omit the reference to the 
algebra and write simply a (A) for the spectrum. For self-adjoint operators on 7i we 
denote the extremal spectral values of cr(A) by A_ = mma(A) and A + = maxa(A), 
so that a (A) C [A_, A+]. 

We will analyze the finite section method for multidimensional index sets of the 
form Z d . To that end we define the projection P n in dimension d > 1. We set C n = 
[—n,n] d H Z d , the integer vectors in the cube of length 2n centered at the origin. 
Then the projection P n is defined by (P n y)(k) = X[~n,n] d {k)y{k) = Xc„{k)y{k) for 
k G Z d . The range of P n is a subspace of £ 2 (Z d ) of dimension (2n + l) d and will be 
identified with C5 2n+1 ) . The finite section is then defined to be A n = P n AP n . By 
definition, A n is a (finite rank) operator acting on £ 2 (Z d ), but we often interpret 
A n as a finite (2n+ l) d x (2n+ l) d -matrix acting on C5 2n+1 ) . In particular, by A^ 1 
we understand the inverse of this finite matrix, but clearly A n cannot be invertible 
on i 2 (Z d ). 

We mention that our results could also be formulated with respect to other index 
sets. 

Theorem 1. If A is a positive and (boundedly) invertible operator on t 2 {Z d ), then 
x n converges to x in £ 2 (Z d ). 

Proof. Step 1. Since by hypothesis, a (A) C [A_,A + ] C (0,oo), we have 
A_||P n 6|| 2 < (AP n b,P n b) = (A n b,b) < A + ||P„6|| 2 . 
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Consequently on the invariant subspace P n £ 2 (Z d ) ~ £( 2n + 1 ) d 

<r{A n ) C [A_,A+] 

independent of n. In particular, each A n is invertible on C^ 2n+1 ^ and 
(2-1) supll^U^Al^ll^lU- 
Step 2. Define an extension of A n by 

(2.2) A n = A n + X + (I-P n ). 

Then a(A n ) C [A_, A + ], and all matrices A n are invertible on £ 2 (Z d ). Furthermore, 

A n = A n + A + (i" — P n ) and A n converges to A in the strong operator topology. 
Step 3. (Lemma of Kantorovich). Since 

WA^b-A^bh = \\A~\a- A n )A' l b\\ 2 

(2.3) < sup IK^IU \\(A-A n )A- l b\\ 2 , 

n 

the strong convergence A n ^ A implies that A n converges strongly to A -1 . 
Step 4. Recall A n x n = b n and Ax = b. Then 

\\x-x n \\ 2 = \\A- 1 b-A- 1 b n \\ 2 = \\A- 1 b-A- 1 P n b\\ 2 

(2.4) < \\(A- l -A n ~ 1 )b\\ 2 + \\A n ~\b-P n b)\\2 = I + II- 
The first term goes to zero by Step 3, and the second term is estimated by 

II < sup IK^IU \\b - P n b\\ 2 < Allfe - P n b\\ 2 

n 

and also goes to zero. □ 

The above theorem uses the £ 2 (Z d )-norm, so this is the realm of C*-algebra 
techniques, cf. the work of Bottcher, Silbermann, et al. [4,17]. 

Several questions arise naturally in the context of the finite section method: 

1. Does the finite section method also converge in other norms, e.g., in weighted 
f-norms? 

2. Can we derive quantitative estimates? If the finite section method works, how 
fast does x n converge to xl What conditions on the matrix A and the input vector 
b are required to quantify the rate of convergence x n — > xl 

3. What conditions and modifications are required (if any) to make the finite 
section method work for matrices that are not hermitian? 

For an answer of the first question, we make the following observation: The 
simple argument above extends almost word by word, provided we can show the 
following properties: 

(1) Both A and A' 1 are bounded on t v m , 

(2) sup„ \\A n \\e m ^e m is finite, and 

(3) the finite sequences are dense in l p m . 
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The answers to the other two questions also revolve around the above observation as 
well as on properties of certain involutive Banach algebras, which will be introduced 
in the next section. 

3. A Class of Banach Algebras of Matrices 

To understand the asymptotic behavior of the finite section method on Banach 
spaces, we need to resort to Banach algebra methods. We first consider some typ- 
ical matrix norms that express various forms of off-diagonal decay. Our approach 
is partly motivated by some forms of off-diagonal decay that is observed in vari- 
ous applications, such as signal and image processing, digital communication, and 
quantum physics. A different way of describing off-diagonal decay of matrices (and 
operators) is given by the notion of band-dominated operators [22] . 

Weights. Off-diagonal decay is quantified by means of weight functions. A non- 
negative function v on Z d is called an admissible weight if it satisfies the following 
properties: 

(i) v is even and normalized such that v (0) = 1. 

(ii) v is submultiplicative, i.e., v(k + 1) < v(k)v(l) for all k, I G Z d . 

The assumption that v is even assures that the corresponding Banach algebra is 
closed under taking the adjoint A*. The weight v is said to satisfy the Gelfand- 
Raikov-Shilov (GRS) condition [10], if 

(3.1) lim v(nk)n = 1 for all k G Z d . 

n— >oo 

This property is crucial for the inverse-closedness of Banach algebras, see Theorem|3] 
below. The standard weight functions on Z d are of the form 

v(x) = e ad{x)b (1 + d(x)) s , 

where d(x) is a norm on R d . Such a weight is submultiplicative, when a, s > and 
< b < 1; v satisfies the GRS-condition, if and only if < b < 1. 
Consider the following conditions on matrices. 

1. The Jaffard class is defined by polynomial decay off the diagonal. Let A s be 
the class of matrices A = (a/a), k, I G Z d , such that 

(3.2) \a M \ <C(l + \k-l\)- s Vk,leZ d 

with norm \\A\\ As = sup fcjieZd \a k i\(l + \k- l\) s . 

2. More general off-diagonal decay. Let v be an admissible weight on Z d that 
satisfies the following additional conditions: v _1 G £ 1 (Z d ) and v^ 1 * v^ 1 < Cv^ 1 (v 
is called subconvolutive) . We define the Banach space A v by the norm 

(3.3) \\A\\ao = sup |o w |t;(A;-0> 

3. Schur-type conditions. Let v be an admissible weight. The class A\ consists 
of all matrices A = (aki)k,iez d such that 

(3.4) sup |afcz| v(k — I) < oo and sup |a&| v(k — I) < oo 
keZd iei, d l&d k& d 
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with norm 



4. The Gohberg-Baskakov-Sjostrand class. For any admissible weight v we define 
the class C v as the space of all matrices A = {dki)k,i& d such that the norm 



5. A further generalization is due to Sun [29]. Roughly speaking, Sun's class 
amounts to an interpolation between C v and A v or between Al and A v . Our results 
also hold for Sun's class, but to avoid a jungle of indices, we stick to the simple 
classes defined above and leave the reformulation of our results in Sun's case to the 
reader. 

These Banach spaces of matrices have the following elementary properties. 

Lemma 2. Let v be an admissible weight and A be one of the algebras A s for 
s > d, A v ,Al,C v . Then A has the following properties: 

(a) Both Al and C v are involutive Banach algebras (i.e., B* -algebras) with the 
norms defined in (|3.4|) and ()3.5|) . A v and A s ,s > d can be equipped with an 
equivalent norm so that they become involutive Banach algebras. 

(b) If A G A, then A is bounded on £ 2 (Z d ). 

(c) If A G A and \b M \ < \a kl \ for all k,l G Z d , then B G A and \\B\\ A < \\A\\ A . 
(A is a solid algebra). 

Proof. Properties (a) and (c) are easy and follow directly from the definition of 
the matrix norms. The statements about A s and A v are proven in [16]. (b) is a 
consequence of Schur's test. □ 

Next we study the spectrum of matrices belonging to one of these Banach alge- 
bras. 

Definition 1. We say that A is inverse-closed in £>(£ 2 (Z d )), if for every A G A 
that is invertible on £ 2 (Z d ) we have that A^ 1 G A. 

Our next theorem states that the matrix algebras introduced above are inverse- 
closed as long as v satisfies the GRS-condition. The precise formulation is slightly 
more complicated, because we need to be a bit pedantic about the weights. 

Theorem 3 (Inverse-closedness) . Let v be an admissible weight that satisfies the 
GRS-condition, i.e., lim^oo v (nk) 1 ^ 71 = 1 for all k G Z d . 

(a) Assume that v~ l G £ 1 (Z a! ) and v~ l * v^ 1 < Cv^ 1 , then A v is inverse-closed 
in B(£ 2 (7j d )) . In particular A s for s > d possesses this property. 

(b) Ifv(k) > C(l + \k\) 5 for some 5 > 0, then Al is inverse- closed in B(£ 2 (Z d )). 

(c) C v is inverse-closed in B(£ 2 (Z d )) for arbitrary admissible weights with the 
GRS-property . 



(3.5) 




(3.6) 
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Remark 3.1. The inverse-closedness is the key property and lies rather deep. While 
for C*- (sub) algebras this property is inherent, for Banach algebras it is always hard 
to prove. Inverse-closedness for A s is due to Jaffard [19] and Baskakov [2,3], a 
simple proof is given in [29]. For A v it was proved by Baskakov [3] and reproved in 
a different way in [16]. The result for C v with v = 1 is due to Gohberg-Kasshoek- 
Wordeman [13] and was rediscovered by Sjostrand [24], the case of arbitrary weights 
is due to Baskakov [3], the algebra A\ was treated by one of us with Leinert [15]. 
More general conditions were announced by Sun [29]. 

The following properties are well-known consequences of inverse-closedness. 

Corollary 4 (Spectral invariance). Let A be one of the algebras A s , A v , A\, or 

C v and assume that v satisfies the conditions of Theorem^ Then 

(a) o-^(A) = o~(A) (the spectrum in the algebra A coincides with the spectrum of A 

as an operator on £ 2 (Z d ) ) 

(b) If A is bounded on l v m for all A G A, then the operator norm satisfies 

(3.8) Pllc-C < C\ a \a for all A e A, 
and 

a e?n (A) C a(A) 

(the spectrum is almost independent of the space A acts on). 

Remark 3.2. Statement (a) is equivalent to inverse-closedness, the norm estimate 
in (b) follows from the closed graph theorem, the inclusion of the spectra is an 
immediate consequence of Theorem 

Let us emphasize that in our analysis of the finite section method we only need 
that the algebra A acts boundedly on H? m . In order to understand how the weight 
m depends on the submultiplicative weight used to parametrize the off-diagonal 
decay, let us briefly discuss some sufficient conditions for the bounded action of A 
on l v m . The weights m satisfy slightly different conditions. Let v be an admissible 
weight. The class of v -mo derate weights is 

(3.9) M v = (m > : sup "^^ < Cv(l), V/ G Z d \. 

I k&d m{k) J 

For example, if a, s 6 R are arbitrary, then m(x) = e ad ^ '(1 + d(x)) s is e^ d ^ b (l + 
<i(x))' s '-moderate. 

The explicit examples of Banach algebras discussed above all act on the entire 
range of i v m for 1 < p < oo and a family of moderate weights associated to v. The 
following lemma provides some explicit sufficient conditions on m for A v , A\ or C v 
to act boundedly on l v m . 

Lemma 5. Let v be an admissible weight. 

(a) If A G A\, then A is bounded simultaneously on all £^(Z d ) for 1 < p < oo 
and m G M. v . 

(b) If A G A v and v (k) := v(k)/(l + \k\) s is submultiplicative for some s > d, 
then A is bounded simultaneously on all t v m {7L d } for 1 < p < oo and m G M. Vo . 
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(c) If A e A v , then A is bounded on £™(Z d ). 

(d) If A e C v , then A is bounded on all £^ n (Z d ) for 1 < p < oo and m e 

Proof. For completeness we sketch the easy proof. 

(a) First, let p = 1, c e ^(Z ') and A G ^4*. Then, since m(fc) < Cw(A; - Z)m(Z), 
we obtain 

fcez d 2ez d feez d «ez d 

«GZ d te& k&L d 

Next, let p = oo and cef™ Then, as before 

\\Ac\\ e? o = sup a ki c i m (k) < C sup ^ \ a ki\ \Q\v(k - l)m(l) 

< C[ sup \ci\m(l) ) sup y2 \a k i\v(k - I) = C\\A\\ A i\\c\ 



The boundedness on £^(Z d ) for 1 < p < oo now follows by interpolation. 

(b) and (d) follow from the easy embeddings A v <^-> A\ , C v C .4* and from (a). 

(c) uses the subconvolutivity of i> . Let A e .Au and c e £^°(Z d ). Then, |a^| < 
ll^llx^C^ ~~ _1 an( i l c zl — ll c lks°' l '(0 _1 - Consequently, 

||Ac||<oo = SUp y^OfcjCj 



fcez d 



«ez d 



< PIUIMk- sup ^ x, 1 JSZm v ( k ) < C\\A\\ Av \\c\\ e ~, 

because (t> _1 < Cv(k)^ 1 . □ 

The matrices of the Banach algebras introduced above can be considered as ap- 
proximate banded matrices. This might suggest that it would be sufficient to set 
those entries smaller than some threshold to zero and simply work with banded 
matrices, which are a special case of sparse matrices. At first sight this may seem 
appealing, since invertible banded matrices have inverses with exponentially fast 
off-diagonal decay [9]. However there is an important difference. In many applica- 
tions, cf. [19,25,26,28] the matrix entries do decay off the diagonal, but thresholding 
would still leave us with banded matrices with the number of non-zero diagonals 
easily in the order of several dozens. The theoretical prediction for the decay of 
the inverse of such banded matrices is so slow that it is meaningless for practical 
purposes. The reason is that by resorting to banded matrices we have neglected the 
decay of the entries above the chosen threshold. Thus banded matrices are simply 
not the most suitable model to capture the decay behavior of those matrices and 
their inverses. 
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4. Finite Sections in Matrix Algebras 

We first study the finite sections of matrices belonging to an inverse-closed matrix 
algebra and give a new characterization of inverse-closedness by means of finite 
sections. This is a necessary step in the qualitative analysis of the convergence 
properties of the finite section method on weighted £ p -spaces, but should be of 
independent interest in the study of Banach algebras. 

Let A FS be the set of all finite sections of matrices in A, formally 

(4.1) A FS = {A= (A n ) neN : A n = P n AP n for some A G A} . 

Although A FS is not an algebra anymore, we may define a notion of inverse- 
closedness. 

Definition 2. We say that A FS is inverse-closed if for every sequence {y4 n }„ e N of 
invertible finite sections such that ||Ai|U — C an d ll^n^Up < C, we have that 
H^-n^U — C"i f° r some constants C and C that do not depend oniiGl 

The comparison of inverse-closedness of A and of A F indicates that the transi- 
tion from the infinite-dimensional setting A to the finite-dimensional case A F is 
done by replacing the hidden word "bounded" by "bounded uniformly in dimen- 
sion" . 

The definition of A FS suggests as a next step to consider sequences of arbitrary 
finite square matrices instead of finite sections. To define a norm that is related to 
the ,4-norm, we must assume that the norm || - \\ A can be applied to arbitrary finite 
matrices by defining an appropriate embedding of (C 2n+1 ) d into £ 2 (Z d ). For j G Z d 
and C n = {—n, . . . , n} d C Z d we define an extension of the (2n + l) d x (2n + l) d - 
matrix B to in infinite matrix B J of B, say B J = {B J ) k ^^d such that 



(4-2) {B J ) j+w 




for k,l G C n 
otherwise. 



Then we define the norm of B by 

(4-3) ||B|U= ||B J |U. 

In the big picture, this definition makes sense only when the norm does not depend 
on the embedding cube J = j + C n . This requires an additional property of A. 

Let T h l G Z d , denote the translation operator Tif(k) = f(k — I) acting on 
/ G £ 2 (Z d ). We say that the norm of A is translation-invariant if 

(4.4) ||T_,AT,|U= \\A\\ A VAeA,leZ d . 

Clearly, if the norm of A is translation-invariant, then ||S J m does not depend on 
the cube J = j + C n , and we can apply || • ||^ to finite matrices. From now on, let 
us assume that || • ||^ is translation- invariant. 

Similarly to (|4.1|) and analogous to [17, Section 1.2.2] we introduce the set A F 

by 

A F = {B = (B n ) nm : B n is a (2n + l) d x (2n + l) d matrix and sup ||5 n |U < °°} > 

rteN 
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and we endow A F with the norm 

(4.5) \\B\\ a f = sup \\B n \\ A . 

neN 

If the norm of A is translation-invariant, then A F is a well-defined object. We note 
that A F is a Banach algebra contained in B := 0^L X B((C 2n+1 ) d ). 

Definition 3. We say that A F is inverse-closed if for every sequence {B m } m€Z d of 
finite invertible matrices such that ||-B m ||.4 < C and H-B^llop < C, we have that 
ll-B^m < C , for some constants C and C that do not depend on m G Z d . 

In view of Definition ^ this amounts to saying that the algebra A F is inverse- 
closed in ©~ =1 £((C 2n+1 ) d ). 

The inverse-closedness of A, A FS , and A F depends on the original algebra A. 
For the study of the relations between them we introduce some further natural 
conditions. 

(CI) Weak solidity: For every A G A there is a constant C such that ||^4 n |U — 

C\\A\\ A for all neN. 
(C2) Weak inverse-closedness: For every A G A that is invertible on £ 2 (Z d ) the 

condition sup neN 1 1 ^4. ~ 1 1 1 < C implies that A^ 1 G A. 

Our third condition concerns an infinite matrix B block that is built from blocks of 
finite square matrices {B m } meN by stacking them "along the diagonal". For this 
we choose a sequence j m G Z d such that sequence of cubes J m = j m + C m C 7L d is 
disjoint. Now set 

(4.6) B block = B J ™ , 

mgN 

where B^ 1 is the extension of B m given in (J4.2IK 
(C3) Block norm equivalence: There exist constants C, C > 0, such that for 
every B G A F 

(4.7) C\\B block \\ A < sup \\B m \\ A < C'\\B block \\ A 

mgN 

Remark 4.1. We want to point out that in all settings A, A FS , and A F it suffices to 
show the inverse-closedness property for positive matrices. Therefore, if necessary, 
one could restrict conditions (C1)-(C3) to such matrices. We note that the upper 
bound in (J4.7)) follows already from condition (CI). 

These three conditions are sufficient to show that the concepts of inverse-closedness 
in A, A FS and A F are equivalent. 

Theorem 6. Let A C B(£ 2 (7j d )) be a Banach algebra such that the norm of A is 
translation-invariant and A satisfies conditions ( CI )-( C3). Then the following are 
equivalent: 

a) A is inverse-closed in B(£ 2 CZ d )). 

b) A FS is inverse-closed. 

c) A F is inverse- closed in B := 0^L X B((C 2n+1 ) d ). 
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Proof, c) =>- b) This implication is clear, because each sequence of finite sections 
A n belongs to A F by (CI). 

b) =>■ a) Since a matrix A is invertible, if and only if the matrices A* A and AA* 
are invertible, we may assume without loss of generality that A is positive and 
invertible on £ 2 (Z d ). So assume that A G A is positive and invertible on £ 2 (Z d ). 
We want to show that, if A FS is inverse-closed, then A^ 1 G ^4. 

By condition (CI), we have that ||A n m < C for all n G N and some constant C 
independent of n G N. Moreover, (|2.1j) implies that ||A n ; 1 || p < ||A _1 || p. Since by 
assumption A FS is inverse-closed, we obtain that H^^IU < C for all n G N and 
some C > 0. By condition (C2) we obtain A^ 1 G A. 

a) => c) We argue by contradiction and show that if (c) fails, then so does (a). 
Assume that A F is not inverse-closed. This means that there is a sequence of finite 
invertible matrices {B m } m€N such that ||-B m m and H-B^Hop are uniformly bounded 
in m G N, but sup mgN H-B^IU = °°- 

Consider the corresponding matrix B block as given in (|4.6|) . Then its inverse 
(^B block )~ 1 is a block matrix that corresponds to the sequence {5~ 1 } meZ d. There- 
fore, \\{B block )- l \\ op = sup mgN ll-B- 1 !!^ < oo, so B block is invertible on £ 2 {Z d ). Con- 
dition (C3) implies that \\B block \\ A < C"sup meZd \\B m \\ A < oo, so B block G A. The 
same condition guarantees that \\(B > C sup mgZ d H-B^IU = 00 • Thus 

[B hlock )- x £ A and A cannot be inverse-closed in B(£ 2 {Z d )). □ 

Next we apply Theorem El to the matrix algebras A v , Al, and C v introduced 
in Section 2. As we mentioned in Lemma |2l all these algebras are contained in 
B(£ 2 (Z d )). Moreover, the norms associated to these algebras are translation in- 
variant. Indeed, to check that (|4.4j) holds, we denote the standard basis of £ 2 (Z d ) 
by {zk}k& d and observe that (T^ATje h e k ) = (AT^T^) = (Aei +j ,e k+j ) for all 
j G Z d . Therefore, {T_jATj)ki = (A)k+j,i+j and since the norms of A v , Al and C v 
use only the difference of k and /, they are translation invariant. 

By Lemma |2fc) each of these algebras is solid, so condition (CI) holds for all 
of them. Condition (C3) is more problematic. Since the norm of a matrix in A v 
and Al is defined in terms of its rows and columns, it follows that \\B \\a = 
sup mgZ d ||_B m ||_4. So property (C3) holds for for A v and Al- Condition (C3) fails, 
however, for C v . 

It remains to consider the weak inverse-closedness (C2). 

Proposition 7. Condition (C2) holds for each of the algebras A v and Al- 

Proof. Assume that A G A is invertible on £ 2 (Z d ) and that sup ngN HA" 1 ]!^ = C < 

oo. Recall that A n = P n AP n + X+(I — P n ) is the extension of A n defined in f!2.2|) . 
Clearly, our assumption that HA" 1 ^ is uniformly bounded, implies immediately 

that \\A n H.4 is uniformly bounded as well. 

Since both A and A -1 are bounded on £ 2 (Z d ), A n converges strongly to A -1 in 

£ 2 (Z d ), as we have seen in the proof of Theorem^ Step 3. Therefore, (A n ei, e&) 
converges to (A~ 1 ei,ek) for all vectors of the standard basis of £ 2 {7L d \ 
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Case I: A = A v . Since \\A n \\ A is uniformly bounded, we have that 
II A» 1 |U= sup \(A n l e h e k )\v{k - I) < C. 

k,l& d 



Thus, we obtain that, for every k, I G Z 



\(A l e u e k )\v{k-l) = lim |(A n e,, e jb )|v(A; - < C, 

n—>oo 

and so A E A with ||A |U < C = sup ngN ||A„ H^. 
Case II: A = A\. We use Fatou's Lemma. We have that 

\\A n \\ A = maxl sup V \{A n e h e k ) \v(k-l), sup V] | (A„ e z , e A )|v (A;-Z) > < C. 

I k&L d iei d ) 



iez d fc kez d 

7d 



Therefore, we obtain, for every k G Z° 

V I (A- 1 ^, e fc ) |u(A; - < lim inf V | (AT^, e fe ) |u(fc - I) < C 

i - — * n^oo * — ' 

and for every / € Z d 

y | (A- 1 ^, e fc ) \v(k -l)< lim inf V | (X~\ u e k ) \v{k -1)<C. 



Taking the supremum over I G Z d (or k respectively), we conclude that A 1 E Al 



-l 



and || A 1 \\ A{: < C = sup neN \\A n \\ A i. □ 

Since we have verified that all assumptions of Theorem El are satisfied for A v and 
Al, we obtain the following result. 

Theorem 8. If A is either A v or Al, then A is inverse- closed, if and only if A FS 
is inverse-closed if and only if A F is inverse-closed. 

Remark 4.2. Theorem |H] holds for Al even if v = 1. In this case, Al is the Schur 
class, i.e., the class of matrices that satisfy the Schur test [18] or, equivalently, the 
class of matrices that are bounded simultaneously on all £ p , 1 < p < oo. It seems 
to be an open problem if this algebra is inverse-closed. Theorem |H] reduces this 
problem to an equivalent (and equally difficult) question about finite-dimensional 
matrices. However, if v satisfies a mild growth condition, see Theorem Efb), then 
Al is inverse-closed. 

We need a few more facts before addressing questions of convergence of the finite 
section method. 

Corollary 9. Let v be an admissible weight satisfying the GRS-condition (J3.1j) . 

(a) If A G Al is positive and invertible on £ 2 (Z d ) 7 then sup ngN HA^H^i < oo. 

(b) Assume in addition that v^ 1 G £ 1 (Z a! ) and v^ 1 *v~ 1 < Cv~ l . If A is positive 
and invertible on £ 2 (Z d ) ; then sup ngN H^^Ha, < °°- 
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Proof. Our assumptions on v imply that A £ {A v , Al} is inverse-closed. By Theo- 
rem El A FS is inverse-closed as well. Therefore, to achieve that sup ngN \\A n ||_4 < oo 
it is enough to assure that ||Ai|U and ||^ 1 ||op are bounded uniformly in n £ N. 
This, however, follows from ||Ai|U < \\A\\^ (solidity) and from Step 1 in the proof 
of TheoremUl where we showed that ||^4~ 1 || 0; p < ||A~ 1 || op . □ 

As the block norm equivalence (C3) fails for the algebra C v , we do not know 
whether (C V ) FS is inverse-closed. However, since C v C Al, we have the following 
result. 

Corollary 10. Let v be an admissible weight satisfying the GRS- condition (13. 1|) 
and v(k) > (7(1 + |A;|) e for some e > 0. If A £ C v is positive and invertible on 
£ 2 (Z d ), then sup ngN H^" 1 ^ < oo. 



5. Convergence in t v m 

After the analysis of finite sections in matrix algebras, we are now in a position 
to show that the finite section method converges in weighted £ p -spaces, whenever 
the matrix is in one of the algebras A s , A v , Al, and C v . 

Theorem 11. Let A be one of the inverse- closed algebras A s , A V ,C V! or Al, where 
the weight satisfies the conditions stated in Theorem^ for each case. Assume that 
A £ A is positive and invertible on £ 2 (Z d ) and acts boundedly on l v m . 

Ij ' b £ i v m and p < oo, then the finite section method converges in the norm of 

ov 

Ij b £ l v m and p = oo, then the finite section method converges in the weak*- 
topology. In particular, x n goes to x entrywise. 

Proof. We expand the model proof of Theorem Q and insert the results about 
Banach algebras obtained in Sections 3 and 4. Recall that A n = P n AP n + A + (7 — P n ) 
is the extension of A n defined in (|2.2|) . Throughout the proof C denotes a constant 
that may change from step to step. 

Step 1 in the proof of Theorem ^ remains unchanged and yields that 

a(A n ) C [A_,A+] 

independent of n and that 

(5.1) sup||,4 n || op < Al 1 = ||v4 _1 || op , 

neN 

(where || • || op is the operator norm on £ 2 (Z d )). Since A is positive and invertible 
on f(Z d ) and A is inverse-closed in B(l 2 {!/•)) by our hypotheses on the weights, 
Theorem El guarantees that A^ 1 £ A as well. By Corollary HJb), the inverse A^ 1 is 
then bounded on t v m . Furthermore, for A £ {A s , A v , Al} by Corollary El we know 
that 

sup \\A n < C sup \\A n m = C < oo . 

neN "' neN 
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For A = C v , Corollary [Tul implies that 

sup \\A n H^-vip < C sup || A n ||_4i = C < oo . 

neN n£N 

Step 2. For p < oo, A n converges to A in the strong operator topology on t p m . 
This follows from the inequality 

\\A - P n AP n \\^ m < ||(J - P„,)A|| C ^ C + ||P n A(/ - P„)|| c _ c , 

and the fact that P n f — > / e £^ is equivalent to the density of the finite sequences 

inC- ^ 

Step 3. (Lemma of Kantorovich) . We know that A n = A^ 1 + A^ 1 (J — P n ). 
Since 

WA^b-A-'bW^ = \\A~\A-A n )A- 1 b\\ em 

(5-2) < suplK^H^elKA-^A-^H^ 

?i 

(5.3) < suplK^miKA-Ai^&ll^ 

the strong convergence A n — A on £^ implies that A„ converges strongly to A" 1 
on £^ for 1 < p < oo. 

Step 4. Recall, that A n x n = b n and Ax = b. Then 

||x - x n \\ e v m = \\A- l b - A- l b n \\e; n = \\A- X b - A^PnbW^ 

(5.4) < ll^-^Vll^ + II^V-P^II^^I + II. 

For 1 < p < oo the first term goes to zero by Step 3, and the second term is 
estimated by 

W b - p nb\\e m < cr*\\b - p n by m , 



(5.5) II < sup H4, || c ^P m \\b - P n b\\& < C-'Wb - P n b\\ t 



n 



and also goes to zero. 

For p = oo we prove weak*-convergence. Assume b e I™ = ^i/ TO J an d V £ ^\/ m - 
Then, 

(x - x n , y) = {A- l b - A?P n b, y) = (b- P n b, A^y) + (P n b, (A~ l - A~ l )y). 

the first term tends to zero, because finite sequences are weak*-dense in £™. The 
second term is majorized by || (A" 1 — A^y^i and converges to zero by 

Step 3, Q. □ 

6. Quantitative Estimates 

In Theorem II II we have investigated the convergence of the finite section method 
in the norm of i v m provided that the input vector b is in l v m . For the quantitative 
analysis, we assume that the input is in l v m and we study the convergence in a 
weaker norm. 
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We first work with the algebra A = A v defined by off-diagonal decay of the 
matrices and a subconvolutive weight satisfying the GRS-condition. Recall that 
C n = {—n, . . . n} d is the cube of integer vectors, so that Y^k^c n ■ ■ ■ becomes a tail 
estimate. 

Theorem 12. Assume that A e A v is invertible and b e £^°(Z d ). Set 

^(n) = (X;«(ib)- a ) 1/a . 

Then the finite section method converges in the l 2 (7j d )-norm with the asymptotic 
estimate for the error 

(6.1) \\x - x n \\ 2 < C(p(n) . 

Proof. We write 

x %n = A b A n b n 

= A- 1 (b-b n ) + A-\A n -A)A- 1 P n b = \ + \\ . 

Estimate of I: 

iiiii 2 < n^iiop n6-6 n || 2 = wa^u^ n 2 ) 1/2 

\k\>n 

< ii^iun&ibCE^H 172 

= II^IU ||6||^^(n) . 

Estimate of II: Set z n = APnA^Pnb = AP n A n 1 P n b, then II = A' 1 (P n - I)z n . 
Using Lemma El Corollary El and the obvious fact that ||P n 6||^ < ||fe||^ (true for 
every solid sequence space) we obtain 

lknllfg° < II^IUJIAj II a v IHU°°j 

or the pointwise estimate 

\{z n )k\<Cv{k)- X 

which is independent of n. So 

|| II || 2 = \\A-\P n -I)z n \\ 2 



2 

|2\V2 



< \\A-'\\ op \\{P n -I)z n \ 

= ii^- 1 iu(Ek^i 2 ) : 

< \\A- l \\ op C{Y,v{k)-*) 112 
= 1 1 A- 1 1| opC(p(n). 



□ 
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n 



-s+d/2 



Remark 6.1. If v(x) = (1 + |x|) s , then <p(n) = (E/^cJ 1 + 1^11)^) ^ 
and we recover the result of [26]. 

Theorem ^] can be generalized to other matrix algebras and sequence spaces. 
For this we note that the input b is in a "small" space, but that we measure the 
rate of convergence in a "large" space. The other item is that A and A n have to 
be invertible on both the small and the large space with uniform bounds. 

The rate of convergence will follow from the following tail estimates for the 
embedding of sequence spaces. 

Lemma 13. Assume that l v m C 
weights. Set r -1 = maxjg -1 



p 1 , 0} and 



(6.2) 

Then, \\b- P n b\\ c < <p(n 

Proof. We write b — P n b ■ 
C„. Then, 116 — P„6|L« = 



for 1 < p, q < oo and m, w are moderate 
w(k) 7 



E 

. k(t.C n 



m(ky 



= (1 — Xc„)b, where \c n is the characteristic function of 
||6m(l - XcJ^h"- If P < Q, then \\c\\ £q < \\c\\& and so 



\b-Pnb\\ { 



< 



bm(l - xc n 



i W 
m 



< 



CP 



I 1 - Xc n 



i W 
m 



I bm I 



(p(n 



(since for r = oo formula (|6.2|) has to be interpreted with the supremum norm). 
If p > q, then r = (g -1 — p^ 1 )^ 1 > 1. Thus, we use Holder's inequality ||6c|| 9 < 



I ell r- and obtain 



\b-P n b\\ C < 



I 1 - Xc„ 



) W 
m 



\\bm\\ip = ip(n) 



□ 



Theorem 14. Let A be one of the inverse- closed algebras A s , A v , A\ orC v . Assume 
that l v m C £^ and t/iat A acts boundedly on both t v m and t q w . If A G A is invertible 
on £ 2 and b G t v m (the "smaller" space), then the finite section method converges in 
t q w (the "larger" space) with the error estimate 

(6-3) \\x-x n \\ tl < C\\b\\ e p m p(n), 

where C = (1 + H^^m) and <p(n) is as in ()6.2|) . 

Proof. As in the proof of Theorem we estimate the error by 

(6.4) \\x - x n \\ c < \\A-\b - b n )\\ez + \\A~\P n - I)z n \\ c , 



where z„ 



AP n A^P n b. 



Since A G A is invertible on £ 2 , by inverse-closedness A 1 £ A and consequently 
A^ 1 is also bounded on t q w . We note that, by Corollaries I§1 and ITU1 we also have 



sup neN \\A n l \\ A < oo. 



For the first term in ()6.4|) we obtain, with Lemma PHH that 
\A-\b-b n )\\ c < ||i4- 1 || 4t ||6-^|| li < WA-^Wl 



\t* m v{n). 
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The second term is estimated by 

\\A-\Pn- I)z n \\ el < \\A- l \\ tl ^ tl \\(P n - I)z n \\ c < || A -1 !!^^ \\z n \\ e p m <f(n). 
Finally, 

Iklk = WAPnA^PnbW^ < \\A\\ Pm ^\\A~ l \\^ f vJ\b\\^ 

and this expression is uniformly bounded by Corollary |§] and El and Corollary 0] 
Thus, we are done. □ 

If 1°^ C £ 2 , then we recover the simpler statement of Theorem [T21 

7. Non-Symmetric Finite Section Method for Non-Symmetric 

Matrices 

In the previous section we derived quantitative estimates for the convergence 
of the finite section method under the assumption that the matrices are positive 
definite. This assumption is crucial. For non-hermitian matrices it is already a 
difficult problem to derive merely qualitative statements about the convergence of 
the finite section method [6,11,12,20]. Indeed, even for very simple non-hermitian 
matrices the finite section method may fail. 

As an example, let us consider the Laurent operator given by the biinfinite 
Toeplitz matrix 



.4 


















1 














c 


1 











c 2 


c 


1 








c 3 


c 2 


c 


1 






Here, we assume |c| < 1 and, as usual in the finite section method literature, the 
box indicates the entry in the zero-zero position. An easy calculation shows that A 
is invertible on £ 2 (Z) and its inverse is the Laurent operator with biinfinite Toeplitz 
matrix 



A' 1 



-c 1 
\*\ 1 

-c 



Let us choose b = eo, i.e, the right hand side is given by the zero-th unit vector. 
The solution to Ax = b is obviously the zero-th column of A' 1 , x = e_i + ceo- 
The finite section method as described in (jl.2|) - (jl.3|) applied to this system fails 
completely, because none of the matrices A n is invertible. In theory the solution 
could be computed by solving the normal equations A* Ax = A*b. Thus, one might 
want to apply the finite section method to the positive definite system A* Ax = A*b 
and invoke the results from the previous sections, since A (and thus A*) belongs 
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to A v ,Al, or C v with weight v(k) = e' fc ' Q ,0 < a < 1. However, the computation 
of P n A*AP n involves the infinite matrices A, A*, which makes this approach not 
feasible for numerical purposes. 

It is easy to see how to alter the finite section method to make it work for this 
particular example. Our goal is more ambitious, and we want to derive a version of 
the finite section method that works for large classes of (algebras of) non-hermitian 
matrices, and not just for some individual cases. We will derive conditions for the 
convergence of the finite section method for non-hermitian matrices in some matrix 
algebras. For this, we consider a slightly generalized version of the finite section 
method. 

Consider the system Ax = b where A is an invertible, but not necessarily her- 
mitian matrix. We set 

(7.1) A r ^ n = P r AP n , and 6 r>n = A* n b, 
and try to solve the system 

(7.2) A r n A r ^ n X r ^ n &r,n 

for properly chosen r and n. Observe that A r<n is a (2r + l) d x (2n + l) d matrix, 
and so A* rn A r ^ n is a (2n + l) d x (2n + l) d -matrix. In general we will need r > n, 
therefore we refer to ()7.1|) - (j7.2J) as non-symmetric finite section method. 

Let us denote B n := P n A*AP n , D Tyn := P n A*P r AP n = A* >n A rjn . Analogously 
to ()2.2|) we define the extensions 

(7.3) B n = B n + \+(I-P n ), D^ n = D r , n + A+(/-P n ), 

where a(A*A) C [A_, A+]. 

Clearly, B n , n G N, is the sequence of finite sections of A* A and B> r n is an 
approximation of B n . We study this approximation for matrices in A v . 

Lemma 15. Assume that A e A v . Then there exists a sequence R(n) G N, such 
that for every r(n) > R(n) 

(7.4) lim \\B n - D r(n) , n \\ Av =0. 

n— >oo 

Ifv(k) = (1 + |A;|) S and A v = A s , then we may choose R{n) = n a for a > an d 
obtain the rate 

\\B n - A^,n|k < C\\A\\ As n a ^ +2s , 

Proof. We define E r ^ n = B n — D rn = B n — D r ^ n . Clearly, E r ^ n is hermit ian and in 
A v , and (E r>n ) k i — for k, I ^ C n . If k, I G C n , then 

{E r , n )kl = {B n )kl — (-Dr,n)fc/ = ^^(^*)fcj a j« — ^^(^*)kj a jl = Q^Oji , 

j& d j'eC r j£C r 

and we obtain the estimate 

\{E r ,n)kl\ < ^ 
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for all entries. If A G A v , then \ajk\ < H^Ha^C/' ~~ Consequently we estimate 
the norm of E r n by 

||-^r,n|U„ = SUp \(E r>n ) M \v(k - I) 



< sup \\a\w <j - krMj - o~M* - 1) 



k,iec n .^ Cr 
Since v(j — l)^ 1 < v(j — k)~ 1 v(k — I), we continue with 

\\Kn\U < \\AWX sup V ^ - k)~ 2 v{k ~ l? ■ 

k > leCn j<£C T 

Clearly, if k, I G C n and j G" C r , then k — / G C 2n and j — G" C r _ n and we arrive 
at the estimate 

hm ; n./j ~ 



(7-5) ||£ nn |k < Pill .sup ^(fc) 2 J2 <fr 

j — n 

As a consequence, we obtain that lim^oo H^nllx = and (|7.4|) is proved. 

If A G A, i.e., v(k) = (l+\k\) s , then sup fceC2 ' n v{kf = 0(n 2s ) and £ Ki)" 2 
0((r — n) d ~~ 2s ). For r(n) = n a , we obtain the explicit estimate 

||S„-D„- n |k <C||Am s n Q ( d - 2s ) +2s , 
which tends to for a > tt^-j. □ 

Theorem 16. Let A G {A,,^} where the weight v satisfies the conditions stated 
in Theorem\^(a). Let Ax = b be given. Assume that b G H? m and that A G A is 
invertible on £ 2 (Z d ) and acts on l v m . 
Consider the finite sections 

(7.6) A* rn A r<n x r 

Then, for every n there exists an R{n) (depending on A_ and v) such that x r r n \ n 
converges to x in the norm of i v m , for every choice r(n) > R(n). 

Proof. We split the error x — x r ^ n into three terms as follows: 

\\x ~ 2r,nlk = \\{A*A)~ 1 A*b - D-lA rin b\\ em 

< \\{A*A)~ 1 A*b - B- l P n A*b\\, Pm + \\B- l P n A% - fl- 1 A,„&|k+ 

(7-7) + WB-'A^b - D-lA r , n b\\ Fm = || I || c + || II || c + || III || c . 

We observe that the vector B~ 1 P n A*b is exactly the result of the finite section 
method applied to the normal equation A* Ax = A*b. Since A* A G A v and A*b G 
l v m . Theorem [TT] is applicable and implies that || I \p m — > for p < oo and I — ► 
weak* for p = oo. 

~-i 

Since A r ^ n = P r AP n and B n P n = B n P n we can estimate the second term by 

|| II Ik = \\B~\p n A*b-P n A*P r b)\\e m 
(7-8) < sup llST'lk-^, ||A*|k||6 - P r 6|k . 
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As in the proof of Theorem ITTl Corollary HI and M imply that sup ngN \\B n \(P m _(P m < 

C sup ngN \\B n \\a v < C < oo. Since the finite sequences are dense in l v m for 
p < oo, ()7.8|) yields || II ||^ — > 0, similarly II — > weak* for p = oo. 
For the third term, we start with the obvious estimate 

II HI ||c = WB-'A^b - D-lA r>n b\\ ePm < \\Bn~ 1 - \\ A *r,nb\U p m ■ 

Here ||y4* n 6||^ = ||P n yl*P r .6|^p i < 11^*11^-,.^ Il^llc ^ s uniformly bounded indepen- 
dent of n and r. 

For the operator norm we use an estimate for inverses in Banach algebras, see 
e.g., [8], and obtain that 

\\B n — D r>n Hc^C — C\\B n — D r>n \\a v 

< g II B n \\av W^n — Dr,n\\A v 



1 ~~ \\B n \\a v \\Bn Dr,n\\Av 



~ 1 __ 

Once again by Corollary El we have sup nGN \\B n \\a v < C < oo, and by Lemma fT5l 

lim^oo \\B n — D r ^ n \\A v = 0. 

Consequently, for any positive sequence e n — > 0, we may choose R(n), such that 
\\B n - D r{n)in \\ Av < e n for r(n) > R(n). 

By combining the estimates for I, II, III, we have thus proved that \\x— a?r(n),n||^ p "~ * 
for every sequence r(n) > R(n) and we are done. □ 

Remark 7.1. If v (k) = (1 + |A;|) S for s > d, then R(n) can be chosen to be n a for 
a > by Lemma HH1 

Remark 7.2. It is well-known that, from a numerical viewpoint, the solution of the 
normal equations should be avoided whenever the condition number of the matrix is 
large. As an alternative to the normal equations one could use matrix factorization 
methods. Since D r ^ n is invertible, the matrix A r>n has full rank (2n + l) d , and one 
could apply a QR-factorization of A r ^ n or some other factorization and compute an 
approximate solution to Ax = b in that way. This idea raises a number of interesting 
questions: For instance, assume we can factorize a matrix A 6 A into A = QR, 
where Q is unitary and R is upper triangular, do the individual components Q 
and R also belong to AP. How about other matrix factorizations such as LU- or 
polar-decomposition? We refer the reader to [27] for answers to these questions. 



We return to the example in the beginning of this section. Clearly, A belongs 
to A v for every weight v(k) = e' fc ' a ,0 < a < 1. Since the entries of A decay 
exponentially off the diagonal, it is not difficult to see that it is sufficient to choose 
r(n) = sn for a sufficiently large s > 1, independently of n. In this particular 
example it would even suffice to set s — n + 1, but as pointed out, our goal was to 
derive a finite section technique that is applicable to large classes of matrices, not 
just to this particular one. 

In light of Theorem ^] it is worthwhile to recall that a necessary and sufficient 
condition for the applicability of the finite section method (jl.2|) - ()1.3j) to Laurent 



QUANTITATIVE ESTIMATES FOR THE FINITE SECTION METHOD 



21 



operators is that the winding number of the invertible Laurent operator is zero, 
cf. [12,17]. For the non-symmetric finite section method the winding number is not 
relevant, the key property is the off-diagonal decay of the matrix. Thus Theorem ITB1 
considerably enlarges the range of applicability of finite section type methods even 
for the classical and thoroughly analyzed cases of Laurent and Toeplitz operators. 



In this section we demonstrate the practical relevance of the theoretical frame- 
work derived in this paper by analyzing a problem arising in digital communication. 
We highlight the details related to the finite section method and refer the reader 
to [21] for a more detailed description of the engineering aspects of the problem. 

In a time-invariant digital communication system one is confronted with a linear 
system of equations Ax = b, where x = {xi\i<zz is a sequence of information symbols 
to be transmitted and b = {bk]k& is the received, discrete signal. We can assume 
that Xk G { — 1, 1}, thus x G The entries of A are of the form 



where h is the channel impulse response, (p is a bandlimited function (the trans- 
mission pulse), T is the transmission period, and R is the receive sampling period. 
We do not go into detail about the particular choice of <p, T, and R. The only facts 
we need are: (i) for properly selected T and R we can choose ip to be a bandlimited 
function in L*(R), where v must satisfy the Beurling-Domar condition, i.e., 



(ii) under certain conditions on h, the matrix A has an inverse for R = T and a 
left- inverse for R < T. 

Furthermore, we note that h is a causal function that decays exponentially in 
time. This implies that A is non-hermitian and that A G A v , the latter follows 
from well-known properties of Beurling convolution algebras [23] and the fact that 
a weight which satisfies (|8.2|) also satisfies the GRS condition (|3.1|) . cf. [14]. 

There are two ways to approach the problem of recovering x from b. In the first 
case we try to recover the entries of x "on the fly", i.e., we solve the truncated sys- 
tem A m ^ n x n = b n . In this case we only assume that b G £°° and the ^-convergence 
estimates of Theorem El apply. In the second case we precompute the inverse 
of A by solving Az = eo where eo is the zeroth unit vector. Due to the specific 
(block)-Toeplitz structure of A, the vector z contains all required information to 
fully determine the inverse of A, which is then used to recover x. In this case we 
apply the non-symmetric finite section method from Section [7| to Az = e®. Since 
A G A v and A*e G t\, quantitative estimates as in Section El apply and we can ap- 
proximate the true solution z with a rate of convergence depending on v. Since in 
this application v can be chosen to be v(x) = eW with a < 1, the (non-symmetric) 
finite section method achieves exponential rate of convergence. 



8. AN EXAMPLE FROM DIGITAL COMMUNICATION 



(8.1) 



&ki = <£>(■ — kR) * h * <p(- — IT) 



(8.2) 




for all 



fc=0 
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